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Robust Dynamical Decoupling Sequences for Individual Nuclear Spin Addressing 

J. CasanovaQ Z.-Y. WangJ^J. F. Haasef] and M. B. Plenicf] 

Institut fur Theoretische Physik, Albert-Einstein-Allee 11, Universitat Ulm, D-89069 Ulm, Germany 

We propose the use of non-equally spaced decoupling pulses for high-resolution selective addressing of nu¬ 
clear spins by a quantum sensor. The analytical model of the basic operating principle is supplemented by 
detailed numerical studies that demonstrate the high degree of selectivity and the robustness against static and 
dynamic control field errors of this scheme. We exemplify our protocol with an NV center-based sensor to 
demonstrate that it enables the identification of individual nuclear spins that form part of a large spin ensemble. 


I. INTRODUCTION 
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The quantum control and detection of individual con¬ 
stituents of nuclear spin ensembles represents a major techno¬ 
logical challenge whose solution would enable the realization 
of robust large scale quantum registers that are required for 
quantum information processing as well as the observation of 
the structure and dynamics of biomolecules. The solution re¬ 
quires both the hardware of a physical sensor that is capable, 
in principle, of detecting the minute magnetic fields emanat¬ 
ing from individual nuclear spins and the software, i.e. the 
control schemes and measurement protocols, that can isolate 
the signal of a target spin from both, the other nuclear spins in 
the ensemble and from environmental noise. 

Recently, the exceptional properties of the nitrogen- 
vacancy (NV) center in diamond have led to rapidly growing 
interest into its use as a quantum device. Microwave radiation 
can manipulate coherently the electron spin of the NV center 
while its initialization and readout can be achieved optically 
mm. These tools enable the control of the hyperfine cou¬ 
pling of the NV center to nearby 13 C nuclei to carry out their 
detection and polarization 00 as well as the realization of 
13 C-NV quantum gates GDIS). Furthermore, NV centers can 
be implanted close to the diamond surface where they can de¬ 
tect the signal of nuclear spins above the surface with single 
spin sensitivity Col suggesting the possibility of examining 
the structure of biomolecules by means of carefully designed 
protocols lU-S]. 

These applications need to achieve simultaneously the de¬ 
coupling from environmental noise while preserving the in¬ 
teraction with the object of interest. This can be achieved by 
dynamical decoupling (DD) protocols that refocus undesired 
couplings by means of control sequences that design filters 
|[T4l - [l8ll transmissive exclusively for very specific interactions 
and frequencies. As small errors in DD protocols can accu¬ 
mulate and exceed the environmental perturbation they are 
designed to protect against, robust schemes as exemplified by 
the CPMG 021 [20), the XY-family of pulses (21] [22) or con¬ 
catenated DD 1231 are highly desirable. Whereas these con¬ 
cepts are well understood for protecting single isolated qubits, 
their extension to multiple and interacting spins to address in- 
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FIG. 1. (Color online) a) and b) Robust composite pulses with adapt¬ 
able time variables x t and y t to increase the spectral resolution. Each 
bar corresponds to a n rotation around the axis defined by tp [denoted 
as (7z%]. The X, a), and Y, b), composite pulses will be connected 
in the robust AXY-n sequence, c) first Fourier component, a half¬ 
period, of the modulation function F(t) for a composite pulse. A 
whole period r corresponds to two successive X and (or) Y com¬ 
posite pulses, d) non-robust pulse unit X where all n rotations are 
applied in the same direction 

dividual spins in an ensemble in the presence of environmental 
noise, remains an outstanding challenge. 

Here we present a method for the tunable control and selec¬ 
tive individual addressing of nuclear spins via the application 
of microwave radiation. Our protocol, denoted adaptive XY- 
n (AXY -n), based on non-equally spaced decoupling pulses 
(241 . generates highly selective filter functions which ensures 
a precise discrimination of the magnetic response of external 
nuclear spins. We demonstrate the robustness of our method 
under realistic conditions including frequency detunings and 
amplitude fluctuations of the microwave driving. In this man¬ 
ner, our protocol represents a realistic procedure to achieve 
high-fidelity quantum gates between the quantum sensor and 
a spin cluster for quantum computing or quantum simulation 
purposes. At the same time, this proposal constitutes a highly 
selective detection method of different nuclear spins in the en¬ 
vironment for structure determination with potential applica¬ 
tions in biology as well as solid state physics. 
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II. BASIC MODEL 

In the presence of a strong magnetic field B z , aligned with 
the NV axis 2, the Hamiltonian of the coupled NV center- 
nuclei system reads (Ti = 1) 

H = DS\ - y e B z S z - y,B z F + S z J] V ?j + H e . (1) 

j j 

Here y e j is the electronic, nuclear, gyromagnetic ratio and H c 
describes the external microwave control fields. The interac¬ 
tion between NV center and nuclei is mediated by the hyper- 
fine vector Aj = \z - 3 ^ft ], fj being the vector con¬ 
necting the NV center and each nucleus in the environment. 
Due to the large zero field splitting D = (27t)2.87 GHz, we 
have eliminated non-secular components in Eq. 0- For the 
sake of clarity we neglect dipolar interactions between nuclei 
in our analytical work, but retain these interactions in our nu¬ 
merics. In the rotating frame of DS 2 - y e B z S z and the control 
field H c we obtain 

H = -Y J <oja> j .r j+ n ± F(t) cr z J] Aj ■ fj. (2) 

j j 

Here, cbj is the unit vector of coj = JjB z z - ^Aj, ojj = \ajj\ 
is the effective Larmor frequency, and cr z = \m s ){m s \ - |0)(0| 
the Pauli operator in the manifold of NV electron spin states 
m s = ±1 and m s — 0. For example, when the m s = 1 electron 
spin state is chosen, a qubit between the m s = 1 and m s = 0 
energy eigenstates is available while the m s - -1 does not 
participate in the dynamics governed by Eq. 0- 

The modulation function F(t) takes the values +1 (-1) de¬ 
pending on whether an even (odd) number of 7r-pulses (in¬ 
stantaneous in our analytical model) have been applied, see 
Fig. |T] The appropriate tuning of F(t ) becomes crucial for ac¬ 
curate spin detection and control. We assume periodic pulse 
sequences, F(t) = F(t + r) so that F(t) = o fk cos (kaJoDt), 
where co DD = 27r/r and the Fourier coefficients are fk = 
2 f Q F(t) cos(ka)DDt)dt. The number of pulses and the spac¬ 
ing between them control F(t) and therefore the values of the 
fk. For example, CPMG sequences with pulse interval r/2 
have only two pulses in a DD period r, resulting in Fourier 
coefficients /o = 0 and fk = 4(kn)~ l sin (y) for k > 0. 

As the interaction between nuclei and sensor is governed 
by the second term of Eq. ([2]), it is convenient to move to a 
rotating frame w.r.t. - 2 j <^j &>j ■ //. We get 

H lnl = m s F(t)f J] /, ■ A/O, (3) 

j 

where Aj{t) = aj cos(c ojt) + Aj x cbj sin(u;/) + Aj • cbj cbj with 
dj = Aj - Aj • cbj cbj. Now we neglect the last two terms 
in Aj(t) because F(t) is symmetric and oscillates with zero 
mean. Additionally we eliminate the counter rotating terms of 
Hint finding 

Hint = aj\_e^r^ DD ) t + cx .]. ( 4 ) 

k j 



FIG. 2. (Color online) Sensor response (transition probability of an 
NV center coupled to a single nuclear spin) for different pulse se¬ 
quences. The total evolution time is around 1.4 ms, B z = 200 G, and 
microwave /r-pulses with a duration of 12.5 ns are applied. A detun¬ 
ing of A = 1MHz and a microwave amplitude mismatch of 5% have 
been included, a) response for a CPMG sequence (600 pulses), b) 
uses the AXY-8 sequence, while in c) each X, or Y is replaced by the 
X composite pulse. The patterns in b) and c) are the result of the ap¬ 
plication of 3040 decoupling pulses with an interpulse spacing such 
that f\ DD = T1 an d f 2 =f 3 =f 4 = o. In b) we denote, with arrows, 
the location of spurious resonances, see the main text. The red line 
crossing the figure denotes the location of the resonance frequency. 


In the following we identify the conditions that allow to tune 
the control field in order to address selectively the ft-th spin. 
When the k DD -th harmonic matches the Larmor frequency, 
kDD^DD = C0 n , we have 


+ — 8~^ X! f kD ° Ij ‘ [ el(C ° J ^ + cx '] 

j+n 

Z ZV- h [e^-^Jk DD )t + c.c.j. 

k±k DD j 

For a sufficiently strong magnetic field coj ~ yjB z » \Aj\ the 
last line in Eq. ([5]) can be ignored if 

\yjB z \ » k DD \dj\. (6) 

Similarly, the second line in Eq. ^ is eliminated if 

I kOj — co n \ ^ \fk DD \Wj\. (7) 

For sequences such as CPMG, XY-type, the Knill DD [25 , 
26], and their combinations (27] [28) where decoupling pulses 
are equally spaced, the coefficient fk DD = 4sm(^ L )/(koD^) 
is not tunable and the strong interactions a j limit the address¬ 
ability because of Eq. 0- Additionally, Eq. 0 limits the 
largest possible k DD . 

Now we will show how our approach can achieve arbitrary 
tuning of the value of the coefficient fk DD to obtain the effective 















3 



FIG. 3. (Color online) Difference between ideal and error-included 
patterns, a) quantity used for comparing the behavior of distinct se¬ 
quences in the presence of errors - the difference between transition 
probabilities, see arrow, averaged over the frequency interval of the 
figure. In b), c), the difference for the sequences used in Figs.[2]b), 
c), is obtained for different values of the detuning and amplitude er¬ 
ror. d) shows the patterns without (solid line) and with (red crosses) 
time-dependent microwave fluctuation in the control field. 


Hamiltonian 


m s -+ 

^int — fkoo^'z^n ‘ 


( 8 ) 


robustness because it is able to cancel first and second order 
error contributions. This allows to extend our protocol to sev¬ 
eral thousands of imperfect decoupling pulses, see Appendix. 
The symmetric construction yields 



^ (-!)/[(_ 1/_ l]sin(2^x ; ) 

j= u 



, (9) 


which vanishes for even k. As an example we choose / 3 = 0 
and /i = f kDD = f\ DD , and a solution is 


*1,2 


1 

— arctan 
2 n 


±(3/100^ - 12)wi + V3 vv2 
V6 a/ w 2 - 96/i DD W\7l ± W 2 V3W2 


( 10 ) 


with the functions w\ = 4 - f\ DD n and W 2 - wi[960 - 
144 /i DD ^-12(/i D0 7T) 2 +(/i D0 7r) 3 ]. Here any f\ DD n e (-8 cos § + 
4, 8 cos | - 4) is possible. Note that for f\ DD ~ 0.16, the 
solution X\ 2 reproduces the 5-pulse repetition unit in ll29l 
(where all n rotations are applied about the same axis) which 
is designed for identifying coupled nuclei in a weak field 
regime and a semiclassical picture which can deviate from a 
full quantum description I30ll3ll . For equally spaced pulses 
AXY-n converges to the KDD^ - KDD^+zi sequences which 
are robust against imperfect pulses [ 26-28). Additionally, we 
can use higher harmonics, k DD > 1, to reduce the pulse num¬ 
ber or to increase the total evolution time at a given magnetic 
field. For /i = 0 and / 3 = fa DD , we obtain 


Xj = - - —- arctan Aq 2 - 1 , 
J 4 2 n y^ J 


( 11 ) 


and therefore individual control of nuclear spins without the with qj = 4 /[ ^5 + nfo DD + (- 1 ) 7 ] and j = 1 , 2 , and / 3dd e 
necessity of using high harmonics, i.e. large values of koD- -). 

' 7T ’ 7T ' ’ 


III. SHAPING THE FILTER FUNCTIONS 

Since fk depends on the overlap between F(t) and each har¬ 
monic component, see Fig.[l]c), we can select sequences that 
tune fk at will. Figs. [T] a) and b) present an example of our 
proposed pulse scheme where each X or Y composite pulse 
is made of 5 7 r-rotations applied along different directions for 
enhanced robustness. When equally spaced, X and Y con¬ 
verge to the Knill composite pulse used in the literature ll25tt . 
Every period r includes two composite pulses X and Y that 
are applied time symmetric around the midpoint of the period 
but are otherwise arbitrary. This guarantees the even charac¬ 
ter of the Fourier expansion of F(t). Accordingly, we have 
5 controllable relative time variables Xj, or yj, each of them 
measuring the relative positions of the first 5 pulses (applied at 
0 < XjT < |). We can use these 5 variables to achieve /o = 0, 
and fkr = 0 with j = 1,2,3 and a tunable coupling constant for 
the ktl l harmonic fk = fk DD • For example, one can adjust the 
interpulse spacing in such a way that /o = /2 = h — fa = 0 
and /i ^ 0, meaning k[ = 2, k[ =3, k[ = 4 and koD = F 

We note that the symmetric construction x 3 = X 4 = 
\ - X 2 , and X 5 = \ - x\ [see Figs. [I] a) and |T] b)] enhances 


IV. ROBUSTNESS AGAINST ERRORS 

In Fig. [2] we compare the transition probability between the 
\x+) = -^=(|0) ± |1)) NV states for different decoupling se¬ 
quences. Possible error sources are included in the control 
Hamiltonian 

H c = Q(t) cos[(tUArv + A)f][5jcCos^(0 + 5^ sin <£>(/)], (12) 

where cojyv is the resonant frequency of the NV transition, 
and the microwave pulses are imperfect, with a detuning error 
A = 1MHz and a 5% of error in the Rabi frequency £l(t). The 
detuning error A takes into account the energy of the nitrogen 
spin intrinsic to the NV center. Because both the magnetic 
field and the hyperfine field Aj for the nitrogen spin are ori¬ 
ented along the NV symmetry axis, flipping of the nitrogen 
spin can be neglected by secular approximation. When the in¬ 
trinsic nitrogen spin is polarized, the detuning A can be made 
small by calibration. If the intrinsic nitrogen spin is not po¬ 
larized, uncertainty of the nitrogen spin state introduce large 
energy shifts on the splitting between the NV states through 
the hyperfine interaction, and the unknown energy shifts cause 
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the detuning error within the range of about 1 MHz for natu¬ 
rally occurring 14 N spin dim 


In Fig. [2] a) the resonance pattern for a CPMG sequence 
with f\ DD = 4 In is plotted. Figs.[2]b) and c) show the result of 
a pulse sequence with fi = f loD = ^ \ and f 2 = f 3 = U = 0- 
While the spectrum in Fig. [2] a) is distorted because of accu¬ 
mulated errors, Fig. [2] b) shows a clear pattern obtained by 
the AXY-8 sequence (AXY-rc for n = 8, i.e. XYXYYXYX), 
where the pulses are repeatedly applied according to the pat¬ 
tern of XY-8. This is a sequence containing 8 composite 
pulses, X or Y, where X and Y are designed according to 
the scheme in Figs. [T] a) and b) respectively. Our simula¬ 
tions extend to 3040 pulses corresponding to a final evolu¬ 
tion time « 1.4 ms. We observe that the resonant peak is 
very sharp which is of considerable benefit for the filtering of 
noisy signals. Note that common to all pulsed schemes are 
spurious resonances due to the finite width of the microwave 
7r-pulses (32). We would like to stress however that the com¬ 
posite pulse structure of the AXY-8 sequence produces very 
narrow spurious peaks which, additionally, are far away from 
the principal resonance even in the presence of noise and im¬ 
perfections (see Fig. [2]b)). Figj2]c) has been computed re¬ 
peating the X sequence, see Fig. [l]d), where all n rotations are 
applied along the same axis. We can observe how the accumu¬ 
lated errors distort the pattern destroying any readable signal. 


To investigate the behavior of each multipulse sequence in 
presence of strong error sources we compare the different pat¬ 
terns that appear when changing the conditions in H c . In Fig. [3] 
a) the cases for A = 0 and perfect microwave pulses (blue 
line), and A = 1 MHz and a 10% of error in the rotating Rabi 
frequency (red line) for the the AXY-8 sequence are plotted. 
In Figs. [3] b) and c) we show the average difference between 
the transition probabilities in the ideal case and in the pres¬ 
ence of error sources in a region close to the resonance posi¬ 
tion. The difference between patterns is small for AXY-8 in 
a detuning range of ±1.5 MHz and ±10% in the microwave 
driving error, see Fig.[3]b). In the case of the combined se¬ 
quence X, Fig. [ 3 ] c), the pulse errors rapidly accumulates. In 
this manner we confirm that the AXY-8 sequence is robust 
against realistic error sources allowing, at the same time, to 
accurately filter specific frequency-components. 


Our AXY-8 is also robust against amplitude fluctuations 
in the driving field. In Fig. [3] d) we compare the re¬ 
sult with (red crosses) and without (solid line) fluctua¬ 
tions in the control field. Fluctuations are modeled by an 
Ornstein-Uhlenbeck processes where the amplitude fluctua¬ 
tion A0(0 evolves according to A £l(t + At) = A£l(t)e~ At/Tmw + 

hq ^ Cmw 2 Tmw (l - g~ 2Af / T mw) (33) . r mw - 1 ms being the correla¬ 
tion time of the microwave noise and c mw « 2d^D 2 r“ 1 w with 
the relative amplitude fluctuation Sq « 7 x 10 -3 , see Ref (231. 
Fig. [3] d) demonstrates that the effect of microwave fluctua¬ 
tions is negligible in the resonance pattern which confirms the 
robustness of our protocol. 
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FIG. 4. a) (Color online) Resonance pattern for a CPMG sequence, 
9 th harmonic used. The position of the vertical lines denote coj/k 
and their height denotes the strength of interaction between the sen¬ 
sor and the nucleus. In b) we can observe the effect of the overlap 
caused by k = 35 (diamonds), 33 (triangles), and 39 (crosses) har¬ 
monics on the spectrum of the k = 37 (circles) harmonic, c) and d) 
show resonance pattern in the AXY-8 sequence for the I st and 3 rd 
harmonics respectively. In each plot, the horizontal axis corresponds 
to the frequency range of the modulation function one should explore 
in order to obtain the nuclear response in the appropriate harmonic. 


V. SINGLE SPIN ADDRESSING 


Finally, we have performed simulations in a diamond con¬ 
taining 736 13 C spins which are generated on random lattice 
sites of the diamond structure according to the natural abun¬ 
dance of 1.1% and B z = 200G. In the simulation we use 
instantaneous pulses and a disjoint cluster expansion includ¬ 
ing dipolar interactions up to a group size of six 13 C (34). 
In Figs. [4] a) and b) the patterns were obtained by CPMG 
sequences with a total evolution time ~ 1.4 ms. Fig. [4] a) 
shows the result for the 9 th harmonic of the 66 -pulse CPMG 
sequence. In this sample, working at higher harmonics causes 
spectral overlap and difficulty in distinguishing signals. This 
can be observed in Fig.|4]b) where the 31 th harmonic (hence 
16 pulses) has been used. Using the AXY-8 sequence, for 
the first [Fig. |4]c)] and third [Fig. [4]d)] harmonic, more spins 
can be resolved by individual resonant peaks. The sequence 
in Fig. |4] c) has been tuned to f\ DD = 37 ^ and fk = 0, k be¬ 
ing 3 or any even number, and it uses 3000 decoupling pulses 
to obtain a total time duration equaling that of Figs. [4] a) and 
b). In Fig.Hd), the third harmonic version with fa DD = 
and fk = Ojc being 1 or any even number, can provide 3-fold 
evolution time (i.e., 4.2 ms) compared to the first harmonic 
version and hence narrower signal peaks. In both cases the 
resolution is increased in a way that permits individual ad¬ 
dressing of a considerable number of spins without the prob- 
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lem of spectral overlaps. 


B. Robust composite pulses 


VI. CONCLUSIONS 


We have presented a protocol that allows for individual ad¬ 
dressing of nuclear spins even in large assemblies and when 
weakly coupled to a sensor. Based on the application of non- 
uniformly spaced decoupling pulses our method obviates the 
need for high harmonics which reduces significantly peak pol¬ 
lution and consequently enhances both resolution and robust¬ 
ness. The method is generic as it permits the high-resolution 
detection of arbitrary magnetic moments by means of a quan¬ 
tum sensor. This scheme promises applications in the control 
of mesoscopic nuclear spin based quantum registers in dia¬ 
mond as well as the precise nuclear positioning in external 
probes including biomolecules. 
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VIII. APPENDIX: ROBUSTNESS OF SYMMETRIC AXY-8 

A. Control errors 

To show the robustness of our protocol for leading and sec¬ 
ond order error contributions on the electron spin rotations it 
is sufficient to consider the pulse Hamiltonian 


2 2 

where A is a detuning, Q(t) = Q during the pulse control and 
Q(t) = 0 between pulses. The phase (p(t) of the Pauli operator 


cr ^ = cos (f) cr x + sin 0 cr y , 


(13) 


1. No pulse delay 


The simplest case is that the pulses are applied subsequently 
without delay. A sequence of composite pulses for the X pulse 
has the propagator (without pulse delays) X= ReRqReRqRe. 


Using Eq. (I4j ) and keeping up to the leading order (first order) 
of errors we find X = Oq + 0{tj 2 ), where 


O 0 = icr x e 


(15) 


is the perfect rotation for the X composite pulse. Therefore 
the X composite pulse cancels the leading order error. 


2. Including pulse delay 


We include pulse delays t t between the (i - l)-th and z'-th 
pulses. In the interpulse spacing, the presence of the detuning 
A cannot be neglected because the accumulated effect A t t is 
not small for long pulse delays. With the pulse delays, the 
propagator of the imperfect composite X pulse is 

X=Ue A5 Uoa 4 Ue A3 U 0 a 2 Ue A i , (16) 

where 

At = 2A(fi), 

A 2 = 2A (t2 ~ h), 

A 3 = 2A(/3 + t\ - £ 2 ), 

A 4 = 2A (£4 + t2 — £3 — t\ ), 

A 5 = 2 Afe + h + t\ - U - ^ 2 ), (17) 

and 


U,p X = i(rj 2 a 2 - l)o> - ir]ecr z U Ai + T]6U Ai - rj 2 p 2 U 4 , (18) 


with a 2 = and U 4 = e The imperfect 

composite Y pulse is related with X pulse by 


Y = U z (tt/2) X UUn/2), 


(19) 


with U z (7t/ 2) = ^* 22 . Now Eq. (16) can be expanded in 
different orders in tj. For symmetric sequences we have 


A4 + A5 — A2 + Aj, 

A 4 = A 2 , 

A 5 = Ai. (20) 


switches for different pulses. One can demonstrate that, up to 
second order on small parameters e = ^ and 6 , the latter being 
the microwave amplitude mismatch, the propagator associated 
to the imperfect rotation is 


R(/> = e 


-'(l-^jcr 0 -/ecr, + ^-^jl (14) 


Note that the small quantities 6 and e can be always written 
as S = rjS and e - pe with 77 < 1 , therefore we can make the 
expansion in terms of the rj parameter. 


One can obtain the first order correction 

Oi=rj(i€(r z -S)(U Al -U A3 ). ( 21 ) 

When the pulses are equally spaced we have f 3 = which 
according to Eqs. © implies Ai = A 3 and therefore U Al - 
U A3 = 0 which leads to 0\ = 0. To suppress higher order 
errors we will combine the composite X and Y pulses in XY 
sequences. Additionally we will demonstrate that even when 
the pulses on the inner level of the X or Y sequences are not 
equidistant but distributed according to Eq. ( [2Q| ) the errors up 
to second order in 77 can be cancelled. 
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3. AXY-4 sequence 

When Ai ^ A 3 , that is the pulses are not equally spaced, 
0\ ± 0. Now we show how leading order error cancella¬ 
tion appear because the combination of the composite X and 
Y pulses, Eqs. ( fl 6 | ) and ( [T9] ), in the XY-4 configuration (this 
scheme constitutes the AXY-4 sequence). Up to the first order 
in T] we have 

X = O 0 + Ou ( 22 ) 

Y = U z (n/2)(O 0 + O x )Ul{nl2), (23) 

Now, by simply calculating the first order contribution in 77 to 
AXY-4 sequences (here XYXY, or YXYX) we have that for 
the XYXY sequence that term is 

OoO z 0 (OoO\ + O x O z 0 ) + (O 0 O\ + OiOfrOo&o = 0, (24) 

while for the YXYX reads 

O z q Oq(O z 0 O\ + O\O 0 ) + (O z 0 O x + O\O 0 )O z 0 O 0 = 0, (25) 

where O z = U z ( 7 r/ 2 )OjU z ( 7 r/ 2 y. This means that terms going 
with 77 are cancelled because of the XYXY, or YXYX, con¬ 
struction. However, one can demonstrate that this cancellation 
does not require the application of the conditions in Eq. ( [20] ). 

4. AXY-8 sequence 

The symmetry imposed by our protocol, i.e. the conditions 
in Eq. ( [20| ), allows us to eliminate the second-order error con¬ 
tribution when we combine our composite pulse with the XY- 
8 sequence, namely, the AXY -8 protocol is applied. Note that 
to cancel the second-order error symmetric composite pulses 
are required. Non-symmetric composite pulse only eliminate 
the the leading order error because of the reasons explained in 


the previous subsection. To make the argument complete let 
us additionally introduce second order error contributions in 
the discusion. This reads 

0 2 = -irfie 1 2 3 4 5 6 + ( 5 2 )(^ 0 -| 0 - j , o -| + cr y ) + 77 2 fi 2 (U Al - U A} ) 
+iT] 2 (6 2 - e 2 ){U A , +2 CT’L(T x (T’i + (J~e0~ X (T | U\ 2+ 2 ) 
+2i] 2 (ed)(U A3+2 crEO- x cr z o--L + cr icr z a x cr* t/ A3+2 ) 
Wif + S 2 )(U Ai _ 3 ct» + cr| £/ Al _ 3 ). (26) 

Therefore, up to the second order one can write 

X = O 0 + Oi + O 2 , (27) 

Y = U z (n/2)(O 0 + O l+ 0 2 )Ul(n/ 2). (28) 

Now under the symmetry conditions in Eq. ( [20] ) we have up 
to the second order of errors XYXYYXYX = (/ 0 +/ 2 X/ 0 +/ 2 ), 
where 

fo = OoO z 0 OoO z 0 , 
fo = OqOoOqOo, 

h = OqO\{OoO\ + O x O\ + 0 2 0\) 

+ (O 0 O z 2 + O x O\ + O 2 ^ 0 )O 0 ^ 0 
+ (O 0 O\ + OxO^XOoOl + O x O z 0 ), 

fi = OlO^Ol0 2 + 0 \Ox + O\0o) 

+ (O z 0 o 2 + 0 \Ox + O z 2 O 0 )O z 0 O 0 
+ (O z 0 Ox + OlOoXO^Ox + O\O 0 ). 

(29) 

It can be calculated that the errors up to the second order fof 2 + 
/ 2/0 = 0 is eliminated by the symmetric contruction of the 
AXY-8 sequences. That is, we have the robust sequence 

AXY-8 = XYXYYXYX = I + 0(r (30) 
where we have used that the error-free control fofo = I. 
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